Topological phases and surface flat bands in superconductors without inversion symmetry 
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We examine different topological phases in three-dimensional non-centrosymmetric superconductors with 
time-reversal symmetry by using three different types of topological invariants. Due to the bulk boundary cor- 
respondence, a non-zero value of any of these topological numbers indicates the appearance of zero-energy 
Andreev surface states. We find that some of these boundary modes in nodal superconducting phases are disper- 
sionless, i.e., they form a topologically protected flat band. The region where the zero-energy flat band appears 
in the surface Brillouin zone is determined by the projection of the nodal lines in the bulk onto the surface. 
These dispersionless Andreev surface bound states have many observable consequences, in particular, a zero- 
bias conductance peak in tunneling measurements. We also find that in the gapless phase there appear Majorana 
surface modes at time-reversal invariant momenta which are protected by a Z2 topological invariant. 
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The hallmark of topological insulators and superconduc- 
tors (SCs) is the existence of topologically protected conduct- 
ing boundary modes. The recent experimental observation of 
these edge and surface states in spin-orbit induced Z2 topolog- 
ical insulators in two and three dimensions (uL respectively, 
has lead to a surge of interest and excitement |2] . An exhaus- 
tive classification of topologically protected boundary modes 
occurring in gapped free fermion systems in terms of symme- 
try and spatial dimension was given in Refs. Interest- 
ingly, this classification scheme, which is known as the "peri- 
odic table" of topological insulators and SCs, predicts a three 
dimensional (3D) topological SC which satisfies time-reversal 
symmetry, but breaks spin-rotation symmetry. Indeed, the B 
phase of 3 He is one example of this so-called "class Dili" 
topological superfluid, whose different topological sectors can 
be distinguished by an integer topological invariant. Recent 
transverse acoustic impedance measurements in 3 He-B con- 
firmed the existence of the predicted surface Majorana bound 
state 01. 

However, finding an electronic analog of the superfluid B 
phase of 3 He remains an outstanding challenge. In this paper 
we argue that some of the 3D non-centrosymmetric SCs might 
be examples of electronic topological SCs in symmetry class 
Dili [7]. We analyze the topological phase diagram of these 
systems and demonstrate quite generally that adjacent to fully 
gapped topological phases there exist non-trivial gapless su- 
perconducting phases with topologically protected nodal lines 
(rings). To characterize these gapless lines we introduce a 
set of topological invariants and show that, due to the bulk- 
boundary correspondence, the presence of topologically sta- 
ble nodal rings implies the appearance of dispersionless zero- 
energy Andreev surface states. These topologically protected 
surface flat bands manifest themselves in scanning tunneling 
spectroscopy (STS) as a zero bias conductance peak, a fea- 
ture which could be used as an experimental signature of the 
topological non-triviality. 

In non-centrosymmetric SCs the absence of inversion in 
the crystal structure generates antisymmetric spin-orbit cou- 
plings (SOC) and leads to a mixing of spin-singlet and spin- 
triplet pairing states. These are the properties that give rise 



to topologically non-trivial quasi-particle band structures in 
these systems. Starting with CePtsSi iJ, a multitude of non- 
centrosymmetric SCs has recently been discovered, including, 
among others, I^Pd^Pts-^B Q9D . 

Model Hamiltonian As a generic phenomenological de- 
scription applicable to any of the aforementioned materi- 
als we employ a single band model with antisymmetric 
SOC and treat superconductivity at the mean field level. 
Thus, let us consider U = ^2 k ^k H ( k )^k with ^ k = 

(cfcf C H, c -fcf c -fc|) T ' w ^ ere °lk * s me e l ectron creation 
operator with spin a and momentum k and the Bogoliubov-de 
Gennes (BdG) Hamiltonian is given by 

H(k) _(Hk) A(k) \ 

The normal state Hamiltonian h(k) describes non-interacting 
electrons in a crystal without inversion center h(k) = Sk&o + 
7fc • <r, where Sk = S-k is the spin-independent part of the 
spectrum, 0x2,3 stand for the three Pauli matrices, and ao de- 
notes the 2 x 2 unit matrix. The second term in h(k) represents 
an antisymmetric SOC with coupling constant 7&. 

Due to the presence of the parity breaking SOC 7& the 
order parameter in Eq. (Q]) is in general an admixture of 
spin-singlet ipk and spin-triplet dk pairing states A(fe) = 
(i^kO'o + dk cr) (i<r 2 ) , where ^k and dk are even and odd 
functions of k, respectively. The direction of the spin-triplet 
component dk is assumed to be parallel to 7^, as for this 
choice the antisymmetric SOC is not destructive for triplet 
pairing floll . Hence, we parametrize the d- vector and the SOC 
as dk = A t Zfc and 7*. = alk, respectively. For the spin- 
singlet component we assume s-wave pairing ipk = A s and 
choose the amplitudes A t?s to be real and positive. 

In order to exemplify the topological properties of 
the BdG Hamiltonian (Q]), we consider a normal state 
tight-binding band structure on the cubic lattice Sk = 
t\ (cosk x + cos k y + cosk z ) — /i, with hopping amplitude 
t\ and chemical potential fi. We will set (t±, /i, a, A t ) = 
(4.0, 4.8, 1.0, 1.0) henceforth. The specific form of the SOC 
7&. depends on the crystal structure lllll . i.e., g~f 9 -ik = Ik, 
where g is any symmetry operation in the point group Q of 



2 




O (b) 



1 2 
' 1 






■J 








v = o : 


: -i 








1 -2 









2.5 



0.5 



1 1.5 



2.5 



FIG. I. (color online). Phase diagram as a function of singlet pair- 
ing amplitude A s and SOC gi [see Eqs. (0 and ©J for the point 
group (a) O and (b) Ca v . The gapped phases are characterized by 
the winding number v with v — (white), z/ = ±1 (dark/light blue 
hatched), v = — 5 (light brown hatched), and z/ = +7 (dark brown 
hatched). Grey shaded and dotted regions are nodal superconducting 
phases with N Cl = ±1 (red/black dotted), Nc 2 = ±1 (light/dark 
grey), Nc 3 = +1 (white dotted), and Nc 4 = +1 (grey). 



the crystal. Having in mind I^Pd^Pta-^B, we assume for 
the pseudo vector 1^ the following form compatible with the 
symmetry requirements of the cubic point group O 




sin k x (cos k y 
sin k y (cos h 
k z (cos k x 




(2) 



with the constant g<i , and where we neglect higher order terms. 
Furthermore, we also consider the point group C± v , relevant 
for CePtsSi, in which case reads 



Ifc — (sin ky&x sin 

+g2 sin k x sin k y sin 



y) 



k z (cos k x 



- cos ky) e z 



(3) 



It is important to note that the quasi-particle band topology 
of H(k), as defined by Eq. (Q]), is mainly determined by the 
momentum dependence of along the Fermi surface sheets. 
Hence, the results we obtain are expected to remain qualita- 
tively unchanged upon inclusion of further-neighbor hopping 
terms in the band structure 

Topological Invariants To characterize the topological 
properties of H(k) we introduce three different topological 
invariants. But before doing so, we observe that H(k) satis- 
fies both time-reversal symmetry (TRS), with T 2 = — 1, and 
particle-hole symmetry (PHS), with C 2 = +1, which are the 
defining symmetry properties of symmetry class Dili in the 
terminology of Ref. |3]. Combining TRS and PHS yields a 
third discrete symmetry, the "chiral" symmetry S = TC, i.e., 
there is a unitary matrix S which anticommutes with H(k). It 
is important to note that while both TRS and PHS relate H(k) 
to H T (—k), S is a symmetry which is satisfied by H(k) at 
any given point k in the Brillouin zone (BZ). 

As shown in Ref. |3] topological sectors in the fully gapped 
phases of H (k) are distinguished by the winding number 
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d 3 k 

24tT 2 * 



Tr [(q-^qXq-^qXq^dpq)] ,(4) 



where the integral is over the 1st BZ and q(k) is the off- 
diagonal block of the flat-band matrix of H (k) 1H2II . 

In the nodal superconducting phases the winding number v 
is no longer quantized. However, we can consider H(k) re- 
stricted to ID loops in reciprocal space and define a topolog- 
ical number in terms of a ID momentum space loop integral 
to characterize the topology of the gapless phases. We ob- 
serve that H(k) confined to a generic momentum space loop 
no longer satisfies TRS nor PHS, but it still obeys chiral sym- 
metry S. Hence, H(k) restricted to a loop in the BZ belongs 
to symmetry class AIII [ 3] and its topological characteristics 
are described by the ID winding number 

N c = ^r£^Tr[g- 1 (/c)V^(/c)], (5) 

where the integral is evaluated along the loop C in the BZ. 
Observe that for any closed loop C that does not intersect with 
gapless regions in the BZ, Nc is quantized to integer values. 
If C is chosen such that it encircles a line node, then Nc de- 
termines the topolo gic al stability (i.e., the topological charge) 
of the gapless line iflUfll. 

Finally, we also consider H(k) restricted to a time-reversal 
invariant (TRI) loop £, which is mapped onto itself under 
k — )■ — k. In that case we obtain a ID Hamiltonian satisfying 
both TRS and PHS (i.e., belonging to symmetry class Dili). 
The topological properties of such a ID system are character- 
ized by the following 7L2 invariant 11211 



W c = II Pf b T ( K )] /Vdet [q(K)], (6) 

K 

where K denotes the two TRI momenta on the loop C and Pf 
is the Pfaffian. Note that Wc is either +1 or -1 for any TRI 
loop that does not cross gapless regions in the BZ. 

Topological Phase diagram Numerical evaluation of the 
topological numbers © and Q yields the topological phase 
diagram of H(k), which is shown in Fig. [T] as a function 
of second order SOC #2 and relative strength of singlet and 
triplet pairing components. Fully gapped phases with dif- 
ferent topological properties (i.e., the phases labeled by v = 
±1, 0, —5, +7) are separated in the phase diagram by regions 
of nodal superconducting phases (grey shaded and dotted ar- 
eas). The fully gapped phases with v = ±1 are electronic 
analogs of 3 He-B. The nodal superconducting phases exhibit 
topologically stable nodal rings, which are centered around 
high symmetry axes of the BZ (see Figs. [2^ and(3^). In or- 
der to determine the topological character of these nodal lines 
(and hence of the corresponding gapless phases) it is suffi- 
cient to consider the topological invariant Nc only for loops 
C that run along high symmetry axes. Thus, for the cubic 
point group O we choose the loops Ci : T —> M —> X T 
and C2 : T — >> M — >> R — ^ T, whereas for the tetragonal point 
group C^ v we consider Cs'.T^Z^R^X^T and 
C 4 :T— ^Z— ^A— ^M— ^T. For the cubic point group we 
find that whenever (Nc x ,Nc 2 ) = (±1,0) there are topolog- 
ically stable nodal rings centered around the (100) axis (and 
symmetry related directions). When ( N^ , Nc 2 ) = (0, ±1) 
the gapless lines are oriented along the (111) axis, whereas 
when (Nc ± , Nc 2 ) = (±1, ±1) the rings are located around the 




FIG. 2. (color online) Nodal rings (a) and (111) surface states (c,d) 
for the point group O with (#2, A s ) = (0.3, 0.5). This parameter 
choice corresponds to the red dotted region in Fig. QJ. (b) Topolog- 
ical invariant iV(m), Eq. 0, as a function of surface momentum 
k\\. Grey and dark blue indicate iV(m) = ±1, while light blue is 
^(iii) = 0. (c) Band structure for a slab with (111) face as a func- 
tion of surface momentum k\ with k\ — 0.757T. (d) Energy disper- 
sion of the lowest lying state with positive energy. The color scale is 
such that black corresponds to zero energy. The states at zero energy 
in (c) and (d) are localized at the surface. The flat bands in (c) and 
(d) are singly degenerate (i.e., one branch per surface), whereas the 
linearly dispersing zero mode at the center of the BZ in (d) is doubly 
degenerate. 

(110) direction. (A similar analysis also holds for the group 

Andreev surface states A non-zero quantized value of any 
of the three topological numbers ©, ([5]) and © implies the 
existence of zero-energy Andreev surface states. First of all, 
in fully gapped phases with topologically non-trivial charac- 
ter there appear linearly dispersing Majorana surface modes 
d, [l5l - [l7ll . In order to understand the appearance of zero- 
energy Andreev surface states in the gapless phases, we now 
make use of the topological invariant Nc with a cleverly cho- 
sen loop C. Let us consider Eq. (Q]) in a slab configuration 
with {Imn) face. In this geometry the Hamiltonian iify mn ) re- 
tains translational invariance along the two independent direc- 
tions parallel to the (Imn) surface. Hence, iify mn ) (fey ) can be 
viewed as a family of ID systems parametrized by the two sur- 
face momenta fey = (k\ , k\ ). Since H^ mn ^ (fey ) obeys chiral 
symmetry (but breaks in general TRS and PHS), its topolog- 
ical properties are given by the ID winding number of class 
AIII 

N (imn)(k\\) = ^~j dk±Tr [q-\k)d±q(k)] , (7) 

where k± is the bulk momentum perpendicular to the surface, 
and d± = d/dk±. Note that Afy mn ) is the same as Nc, 
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FIG. 3. (color online) Same as Fig. [2 but for the point group C± v , 
for a slab with (012) face, and with (g 2 , A s ) = (0.0,0.5). This 
parameter choice corresponds to the white dotted area in Fig. [TJd. 



Eq. (0), with C chosen along k±, following a non-contractible 
cycle of the BZ torus T 3 . 

Now, the key observation is that the above line integral is 
closely related the loop integral Nc, with C = d, that de- 
termines the topological charge of the superconducting nodal 
lines. That is, for those surface momenta fey for which the 
loop along k± in Eq. © passes through just one non-trivial 
nodal ring, Afy mn )(fey) is equal to the topological charge of 
this given nodal ring. Hence, if we plot Afy mn ) (fey ) as a func- 
tion of surface momenta (see Figs. EtOJ)), we find that the 
boundaries separating regions with different winding num- 
ber are identical to the projection of the nodal lines onto the 
(Imn) plane. Furthermore, since a non-zero quantized value 
of Afy mn ) implies the existence of zero energy states at the 
end points of the ID Hamiltonian Hn mn ) (fey ) olflsll. we find 
that there are zero-energy Andreev bound states on the (Imn) 
surface located within the projected nodal rings. This conclu- 
sion is corroborated by numerical computations of the zero- 
energy surface states both for the point group O and C± v (see 
Figs. [Hand [3]) . When two nodal rings overlap in the (Imn) 
projection of the BZ, then the quantized value of Afy mn ) in 
the overlapping region is determined by the additive contribu- 
tion of the topological charges of the two rings. In particular, 
one can have a situation where the two contributions cancel, 
in which case there is no zero-energy surface state within the 
overlapping region. 

Finally, using an analogous argument as in the previous 
paragraph, we can also employ the Z2 number © to deduce 
the presence of zero energy modes at TRI momenta of the 
surface BZ II 1211 . One example of this is the Kramers pair of 
surface zero modes located at the center of the surface BZ in 
Fig. Hi (cf. Refs. EE El). Remarkably, this is a surface Ma- 
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FIG. 4. (color online). Surface and bulk density of states for the 
point group (a) O and (b) Cav . The surfaces are oriented perpendic- 
ular to the (111) and (012) axes, respectively. The employed param- 
eter values are the same as in Figs. [2] and [3] 



jorana mode in a gapless (nodal) superconducting phase IIS 

Experimental signatures One of the most direct signatures 
of the topological aspects of non-centrosymmetric SCs are the 
surface Andreev bound states. These can be probed by angle- 
resolved photoemission measurements, or by STS of the sur- 
face density of states (SDOS). STS has proved to be an ef- 
fective tool to explore surface states of two-dimensional un- 
conventional superconductors, see, e.g., i20l - [24ll . The bulk 
density of states of 3D gapless SCs with nodal lines vanishes 
linearly at zero energy. In contrast, the surface flat bands lead 
to a diverging zero-energy peak in the SDOS (see Fig. [4]). 

The zero-bias peak in the SDOS is strongly dependent on 
the surface orientation. From this dependence it is in prin- 
ciple possible to (partially) map out the location of the topo- 
logically stable nodal lines in the bulk BZ. In addition, one 
can take advantage of the fact that an applied magnetic field 
leads to a splitting of the zero-energy peak. Again, this split- 
ting is strongly dependent on the orientation of the magnetic 
field axis with respect to the nodal lines. Another possibil- 
ity is to use spatially resolved STS to investigate the SDOS 
in the presence of impurities on the surface. It is expected 
that surface impurities will lead to strong spatial modulations 
of the SDOS, which might give some information about the 
topological characteristics of the nodal lines in the bulk. 



Stability of surface modes The zero-energy surface flat 
bands in time-reversal symmetric non-centrosymmetric SCs 
are topologically protected against the opening of a gap and 
are therefore stable against weak symmetry preserving de- 
formations. Conversely, any perturbation that leads to a gap 
opening of the surface states is expected to be accompanied by 
the breaking of the symmetries of the time-reversal symmetric 
SC, i.e., TRS or certain types of translational invariance. One 
possible scenario, for example, is that interactions might lead 
to spontaneous TRS breaking at the boundary of the SC, such 
as to the coexistence of TRS breaking and TRS preserving or- 
der parameters near the surface. This would be observable in 
experiments, for instance, as a splitting of the zero-bias con- 
ductance peak. 

In conclusion, using three different topological invariants, 
we examined the topological properties of general 3D non- 
centrosymmetric superconductors with TRS. We showed that 
in nodal superconducting phases there always appear disper- 
sionless Andreev surface bands. We established a correspon- 
dence between these zero-energy surface flat bands and the 
topologically protected nodal lines in the bulk, thereby reveal- 
ing the topological origin of the surface flat band. In partic- 
ular, we demonstrated that the projection of the nodal lines 
on the surface coincides with the boundary of the surface flat 
band. We emphasize that the presented formalism (or a gen- 
eralization thereof) can be applied to any 3D unconventional 
SC that preserves TRS. One particularly interesting family of 
compounds is Li2PcL c Pt3_ :E B. In these SCs the substitution of 
Pd by Pt seems to be related to the relative strength of singlet 
and triplet pairing states i25ll . Hence, it might be possible to 
observe in I^Pd^Pts-^B the transition between two topolog- 
ically distinct quantum phases as a function of Pt concentra- 
tion. 
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Supplementary Materials 

We first discuss basic symmetry properties of superconduc- 
tors with time-reversal invariance and then go on to derive the 
topological numbers ©, ©, and © from the main text. We 
shall keep the analysis as general as possible, such that it may 
be applied to arbitrary superconducting systems. In SectionlDl 
we will then specialize to the Bogoliubov-de Gennes Hamil- 
tonian (Q]) describing a single-band non-centrosymmetric su- 
perconductor. 

Appendix A: Symmetries of the Bogoliubov-de Gennes 
Hamiltonian 

Let us consider a general time-reversal invariant supercon- 
ductor belo nging t o symmetry class Dili in the terminology 
ofRefs. USE! 

H(k] _fh(k) A(k) \ 
H ^-\A\k) -h T (-k))> (A1) 

with the TV-band normal state Hamiltonian h(k) and the 
superconducting gap matrix A(fe), which obeyes A(fe) = 
—A T (—k) because of Fermi statistics. With TV bands (or- 
bitals) and two spin degrees of freedom, the total dimension 
of the Bogoliubov-de Gennes Hamiltonian at momentum k 
is 47V x 47V. A class Dili superconductor satisfies two in- 
dependent anti-unitary symmetries: time-reversal symmetry 
T = KXJt, with T 2 = —1, and particle-hole symmetry 
C = JCUc, with C 2 = +1. Here, JC stands for the com- 
plex conjugation operator. Time-reversal symmetry constrains 
H(k) as 

U T H*(-k)U ] T = +H(k), (A2) 

with Ut = diag(i^T, Uj) and ut is a 2N x 2N unitary matrix 
that implements time-reversal invariance of the normal state 
Hamiltonian, i.e., UTh*{—k)v} T = h(k) and = — ut- 
Observe, Eq. (IA2b implies urA^(k) = A(k)uj,. Particle- 
hole symmetry acts on the Bogoliubov-de Gennes Hamilto- 
nian H(k) as 

U c H*(-k)tf c = -H(k), (A3) 



where Uc = &i ® l2JV» 01,2,3 stand for the three Pauli matri- 
ces, and J.2JV is the 27V x 2N unit matrix. Combining time- 
reversal and particle-hole symmetry we obtain a third discrete 
symmetry, which is given by 

UlH(k)U s = -ff(fc), (A4) 

with Us = iUrUc- In other words, there is a unitary ma- 
trix Us that anticommutes with H(k) and thereby endows the 
Hamiltonian with a "crural" structure. Namely, in the basis in 
which Us is diagonal, H(k) takes block off-diagonal form 

H(k) = VH(k)V^ = ( Dt ° (fe) D ^) , (A5) 

where the unitary transformation V is given by 

V=^(] 2N +mT V (A6) 
y/2 \l2N -iu T J 

and the block off-diagonal component reads D(k) = h(k) + 
iA(k)uJp. In the off-diagonal basis the unitary matrix Ut is 
given as Ut = VUtV t = g\ <g> ut. Thus, time-reversal 
symmetry acts on D(k) as follows 

u T D T (-k)u ] T = D(k). (A7) 

To compute the %2 invariant, Eq. ©, it is advantageous to 
perform a second basis transformation which brings Ut into 
the simple form Ut = WUtW t = iai <S> t2N- This 
can be achieved with the help of the unitary matrix W = 
diag(l2Ar, ut)- In this new basis the Bogoliubov-de Gennes 
Hamiltonian reads 

H(k) = WH(k)W^ = ^ t ° (fe) D W) , (A8) 

with the block off-diagonal component 

D(k) = D(k)u ] T = h{k)u ] T - zA(fe). (A9) 

We note that time-reversal symmetry operates on D(k) as 
D T (-k) = -D(k). 



Appendix B: Flat Band Hamiltonian and Winding Number 

For the derivation of the topological invariants it is conve- 
nient to adiabatically deform if (fe), Eq. (IA1I) . into a flat band 
Hamiltonian Q(k). The only assumptions that we need for 
computing Q(k) are: (i) the Hamiltonian has a full spectral 
gap and (ii) there is a unitary matrix Us anticommuting with 
H(k). Thus, the following derivation of Q(k) is applicable 
to any chiral symmetric Hamiltonian with a full bulk gap, in 
particular also to the three-dimensional topological supercon- 
ductors in symmetry class AIII, Dili, and CI O, EL H3-H1 . 

In what follows, we work in a basis in which H(k) takes 
block off-diagonal form. The flat band Hamiltonian is defined 
in terms of the projection operator P(k) which projects onto 
filled Bloch eigenstates of if (fc) at a given momentum k. The 
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projector P(k), in turn, is defined in terms of the eigenfunc- 
tions of H(k) 

D(k)\ f X t(k)\ _ f X i(k)\ 
Dt(fc) ){vt(k))- ±Ao( fe) U(fc)J' ( B1) 

where a = 1, . . . , 2N is the combined band and spin index. 
We assume there is a spectral gap around zero energy with 
|A a (fe)| > 0, and for definitiveness we choose A a (fc) > for 
all a. Multiplying equation (IB lb from the left by H(k) yields 

DD^xt(k) = Xlxt, D^D V t = \ 2 a vt- (B2) 

Hence, the eigenfunctions (xt •> vt) can ^ e obtained from the 
eigenvectors of DD^ or D^D 



DD^u a = X 2 a u a , D^Dv a = X 2 a v a . 



(B3) 



The eigenvectors u a , v a are taken to be normalized to one, i.e., 
u\u a = v\v a = 1, for all a (here, the index a is not summed 
over). The eigenstates of D^D follow from the eigenstates of 
DD^ via 

v a =N a D^u a , (B4) 

with the normalization factor Af a . Using Eq. dB3b one can 
check that v a is indeed an eigenvector of D^D, 

D^Dv a = D^D(M a D^u a ) = M a X 2 a D^u a = X 2 a v a , (B5) 

for all a. The normalization factor Af a is given by 

u\Drfu a = X 2 a u\u a = X 2 a N a = (B6) 



for all a. It follows that the eigenfunctions of H (k) are 



1 ( u a 
V2 \±v a J 72 \±&u a /\ a 



■ (B7) 



With this, the projector P(k) onto the filled Bloch states be- 
comes 



(B8) 



1 /^1 2 JV 

2 V 



^ _ - W ^ 

l2Ar7 2 ^ \p a u\ 



Finally, we obtain for the flat band Hamiltonian Q, which is 
defined as Q = 1 4JV - 2P d, 




(B9) 

In other words, the off-diagonal block of Q(k) reads 

«( fe ) = E y4rTMfc)<4(fc)£(fc), (BIO) 

where u a (k) denotes the eigenvectors of DD^ . For a sys- 
tem with completely degenerate bands, A a = A, for all a, the 
above formula simplifies to 

«( fc ) = A^y E^( fe ) w i( fc ) D ( fe ) = ^y D ( fc )- (B11) 



Examples of topological insulators and superconductors with 
completely degenerate bands are the Dirac representatives of 
Ref. H. 

The integer-valued topological invariant characterizing 
topological superconductors is now simply given by the wind- 
ing number of q(k). It can be defined in any odd spatial di- 
mension. In three dimensions we have 



^3 



BZ 



d 3 k 

24tt 2 



s^Tt [(q-%q)(q-' 1 d 1/ q)(q- 1 d p q)] , 



(B12) 



(B13) 



and in one spatial dimension it reads 

v\ = — / dkTv [q^duq] . 

Alternatively, it is also possible to define the winding num- 
ber in terms of the unflattened off-diagonal block D(k) of the 
Hamiltonian. For example, for the winding number in one 
spatial dimension this reads 

vi = ^~ I dkTr [D~ l d k D - {Z^-i^tl 
4ttz J bz 

= — Im / dkTr[d k \nD}. 
27r Jbz 



(B14) 



Appendix C: Z 2 Invariant for Symmetry Class Dili 

In this section we compute the Z 2 topological invariant for 
symmetry class Dili in d = 1 and d = 2 spatial dimensions. It 
is most convenient to perform this derivation in the basis (IA8lh 
in which the 4 A 7 " x 4 A" Bogoliubov-de Gennes Hamiltonian 
takes the form 



H(k) 



D{k) 
Z)t(fc) 



D(k) = -D T (-k). (CI) 



In this representation, the time-reversal symmetry operator is 
given by T = KXJt = /C 202 I27V and the flat band Hamil- 
tonian reads 

0(fe) = ( /(k) qi o ] ) ' q{k) = -« T (- fc )- ( C2 > 

The presence of time-reversal symmetry allows us to define 
the Kane-Mele Z 2 invariant 11001, 



Pf [w(K) 



K \/det [w(K)}' 



(C3) 



with If a time-reversal invariant momentum and Pf the Pfaf- 
fian of an anti-symmetric matrix. Here, w(k) denotes the 
"sewing matrix" 



w ab (k) = (u+(-k)\Tu+(k)), 



(C4) 



where a,b = 1 , . . . , 2N and w„ (fe) is the a-th eigenvector 
of Q(k) with eigenvalue ±1. The Pfaffian is an analog of 
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the determinant that can be defined only for 2n x 2n anti- 
symmetric matrices A. It is given in terms of a sum over all 
elements of the permutation group S2 n 



Pf(A) 



2 n n! 



sgn(o-) JjA^i-i),^*). 



CTES2n 



Due to the block off-diagonal structure of Eq. (IC2t a set of 
eigen Bloch functions of Q(k) can be constructed as 



or, alternatively, as 

l«^(*0>s 



J_ / n 



1 / ±q(k)n a 
V2\ n a 



(C5) 



(C6) 



where n a are 2 AT momentum independent orthonormal vec- 
tors. For simplicity we choose (n a )& = 5 a b. In passing, we 
note that both |i^(/c))n and |^(fc))s are well-defined glob- 
ally over the entire Brillouin zone. To compute the 7L 2 topo- 
logical number we choose the basis \uf(k))^. Combining 
Eqs. CU) and O yields 



Wab(k) 



"5( 



n\q(—k) )i(T2 (8) J-2iv/C 



q T (k)n b 



n b 



q*(k)n b 



rib 



I {nlq T (k)n b -nlq(-k)n b 



(CI) 



In the second last line we used Eq. (IC2K i.e., q(—k) = 
—q T (k). In conclusion, the Z2 topological number in spatial 
dimensions d = 2 and d = 1 is given by 



W 



n 



Pf [<f (*)] 



L K X Vdet[ 9 (Jf)]' 



(C8) 



where denotes the four (two) time-reversal invariant mo- 
menta of the two-dimensional (one-dimensional) Brillouin 
zone. 



Appendix D: Flat Band Hamiltonian for the 
Non-centrosymmetric Superconductor, Eq. (1) 

Let us now apply the formalism developed in the pre- 
ceding sections to the Bogoliubov-de Gennes Hamiltonian 
(Q} from the main text, describing a single band non- 
centrosymmetric superconductor ll38[|39ll . First, we note that 
time-reversal symmetry for H(k), Eq. (Q]), is implemented by 
U T H*(-k)U ] T = +H(k) with U T = a <8> ia 2 . Hence, 
we need to set ut = icr 2 in Eq. (IA2I) . It then follows from 
Eq. dA6b that H(k) can be brought into block off-diagonal 
form by the unitary transformation 



V 



J_ ( 1 2 -cr 2 
^^12 +cr 2 



(Dl) 




FIG. 5. (a) Band structure of Hamiltonian ([T} for the point group O 
in a slab geometry with (111) face as a function of surface momen- 
tum k\ with k\ — 0, i.e., along the dashed line in panel (b). Here, 
we set (02, A s ) = (0.3, 0.5). (b) Brillouin zone of the (111) surface 
with the values of W(m) (K\\ ), Eq. dD12K at the four time reversal 
invariant momenta K\\ G {r, Mi , M2, M3}. 



The transformed Hamiltonian is given by 
with the off-diagonal block 



(D2) 



D(fe) 



A(l% + ill) B k 
= B k a + Al k -a, (D3) 

and where we have introduced the short-hand notation 

A = a + iA t , B k =s k +iA s . (D4) 

Alternatively, we can also choose to work in the basis dA8l) . in 
which case the off-diagonal component reads 



rVM ( A ( l % ~ O 



-B k - All 



(D5) 



For the computation of the flat band Hamiltonian Q{k) it is, 
however, more convenient to use Eq. dD3l) . 

Repeating the steps of section|Bl we calcualte the eigenvec- 
tors u a (k) of 



D k D{ = 



\B k \ 2 - 



{AB* k +B k A*)l k -a, (D6) 



where l k = \l k \. The eigenfunctions u a (k) of D k D k can be 
obtained by diagonalizing l k ■ cr . Hence, when (l k ,l k ) ^ 
(0, 0), we find that the eigenvectors u a (k) are given by 



1/2O) 



^2l k (l k T 1%) V ±l >< ~ L k 



^11 



(D7) 



According to Eq. (IB11K the off-diagonal block of the flat 
band Hamiltonian Q(k) is defined in terms of the eigenvec- 
tors u a (k). Thus, we need to compute 

-^-u a (k)u{(k) (D8) 

a=l,2 Aak 

(Aifc + A 2 fc)cr + (A 2 fc — Aifc)-^ • cr 



2AifcA2fc 
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with the two positive eigenvalues Xik = \B k — Al k \ and 
^2k = \Bk + Alk\. Note that the last term in the second 
line of Eq. (ID8I) contains removable singularities at the points 
&o where l ko = 0. For those points in the Brillouin zone one 
needs to carefully take the limit k —> ko to obtain the correct 
value of Eq. (ID8I) . Finally, by use of Eq. (IB lib together with 
Eqs. (lD3b and (ID8I) we find for the off-diagonal block of the 
flat band Hamiltonian 



q(k) 



1 



{Al k (X 2 k - Aife) + B k (X lk + A 2 fc)} cr 



{^/^(Aifc + \2k) + B k (\2k — Aife)} — • cr 



• (D9) 



Now, for the 7L 2 invariant we need to bring q(k) into the 
basis in which Ut = icr 2 <S> H-2- This is achieved by letting 



q(k) 

Using Eq. (IC8I) we get 



-iq(k)a 2 . 



Pf [i^(ic)] =n 




(D10) 



(Dll) 



where we have made use of the fact that l k is an antisymmetric 
function, i.e., = — 



1. surface state 



As discussed in the main text, the 7L 2 number (IDllb can be 
used to deduce the presence of Andreev surface states at time- 
reversal invariant momenta of the surface BZ. To exemplify 
this, let us consider Hamiltonian (QJ in a slab geometry with 
(Imn) face. At the four time-reversal invariant momenta K\\ 
of the {Imn) surface BZ the 7L 2 invariant is defined by 

■w»,)=n " li " MT(K± - K ' )] ■ <■>*> 



K ± ^det [ia 2 q T (K ± ,K\\)] 



Eq. (ID 121) is quantized to +1 or —1, with W^ mn ^ (K\\) = — 1 
indicating the presence of Kramers degenerate surface modes 
at the surface momentum K\\ (see Fig. [5]). 
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